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Abstract 

Combined  Direct-Adjoint  Approximations  for  Large-Scale  Design-Oriented 
Structural-Acoustics  Finite-Element  Analysis 

Levent  Coskuner 

Chair  of  the  Supervisory  Committee 
Dr.  Eli  Livne 

Department  of  Aeronautics  and  Astronautics 

A  second-order  combined  approximate-direct/approximate-adjoint  method 
for  accurate  and  efficient  reduced  order  calculation  of  frequency  response  for  a 
structural  acoustic  enclosure  modeled  using  finite  elements  is  presented.  A 
design-oriented  approach  was  adopted:  not  only  does  the  method  account  for 
responses  itself  (such  as  displacements,  stresses  and  pressure  levels)  it  also 
accounts  for  sensitivities  of  the  responses  with  respect  to  design  variables. 

Approximate  first  order  results  were  obtained  using  mode-displacement 
order  reduction,  reduced-order  adjoint  solutions  based  on  Ritz  vectors,  and  first 
order  Taylor  series.  A  combined  second-order  approximation  incorporating  those 
first-order  approximations  was  then  created  for  models  with  modified  design 
variables,  and  was  compared  to  full  order  results  of  a  large  finite-element 
problem.  Accuracy  of  the  second  order  method  was  quantified. 
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Introduction 
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The  examination  of  large-scale  structural-acoustic  finite-element  models 
in  the  time  or  frequency  domain  is  an  important  part  in  the  analysis  of  complex 
physical  structures,  especially  in  the  development  of  design  optimization  tools 
(Refs.  4,  7,  9,  12,  14,  15,  17,  25,  27,  29,  32,  36,  41).  Such  finite-element  models 
often  contain  tens  of  thousands  or  even  hundreds  of  thousands  of  degree  of 
freedoms,  translating  to  a  system  of  equations  of  equally  formidable  size.  Solving 
these  systems  of  equations  multiple  times  for  every  frequency  in  the  frequency 
domain  is  expensive  computationally  (Ref.  16).  For  design  optimization  problems 
where  behavior  sensitivities  are  required,  and  the  system  of  equations  have  to  be 
solved  repetitively  to  drive  the  system  towards  its  optimum,  the  cost  of 
computation  becomes  prohibitive  (Refs.  3, 13). 

It  is  therefore  desirable  to  replace  the  detailed  large-scale  models  with 
computationally  efficient  approximations  of  reduced  order,  while  preserving  the 
accuracy  of  the  simulation.  Commonly  used  order  reduction  methods  in 
multidisciplinary  design  optimization  today  include  first-order  Taylor  series  based 
direct  and  reciprocal  models,  Guyan  reduction  (Ref.  11),  Ritz  vectors  (Ref.  38), 
Lanczos  coordinates  (Ref.  28)  and  Modal  displacement  (Refs.  4,  13,  14).  These 
methods  seek  to  create  a  group  of  reduced  basis  vectors  from  a  reduced  set  of 
shape  functions  that  will  reduce  the  order  of  the  problem.  While  accurate  in 
predicting  deformations,  these  methods  are  often  inadequate  to  accurately  predict 
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stress  and  sensitivity  information,  mainly  because  of  the  need  to  capture  local 
effects  in  such  cases.  Furthermore,  for  structural-acoustic  problems,  methods 
involving  modal  approximations  face  an  additional  hurdle  that  in  most  of  the 
separate  structural  and  acoustic  modes  that  are  needed  to  approximate  the 
problem  are  often  far  apart  along  the  frequency  band  (Ref.  24). 

To  improve  accuracy  of  approximations  obtained  using  first  order 
sensitivity  information,  a  method  combining  solutions  of  the  approximate  direct 
and  adjoint  problems  has  been  shown  to  accurately  predict  stresses  for  airframe 
structures  in  the  frequency  and  time  domain,  as  well  as  unsteady  aerodynamic 
forces  on  wings  in  the  frequency  domain  (Refs.  20,  22,  23).  It  was  demonstrated 
that  this  combined  direct  and  adjoint  approximate  approach  has  second  order 
accuracy  (Ref.  34),  and  is  similar  to  the  second-order  accuracy  shown  by  the 
Rayleigh  Quotient  approximation  (Ref.  5). 

The  Second-Order  Approximation  (SOA)  method  presented  in  this  thesis 
combines  two  first-order  approximations  of  significantly  reduced  order  to 
accurately  calculate  response  and  design  sensitivities  of  a  complex  large-scale 
structural-acoustic  problem  with  thousands  of  degrees  of  freedom,  while  reducing 
this  original  large-scale  problem  into  a  much  more  manageable  and 
computationally  efficient  model  of  merely  tens  of  equations.  In  this  work,  the 
first-order  approximations  used  are  the  Modal  Displacement  (MD)  method  and 
the  Ritz  Vector  (RV)  method.  Performances  of  these  first  order  methods  are 
compared  to  that  of  the  SOA  method,  specifically  when  dealing  with  design 
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sensitivities,  for  which  first-order  methods  are  generally  less  accurate.  Also, 
problems  that  arise  when  using  modal  approximation  methods  for  models  where 
material  properties  differ  by  several  orders  of  magnitude  (for  example,  the  density 
of  air  compared  to  that  of  metal)  or  where  modes  of  interest  are  very  far  apart  in 
frequency  range  (for  example,  beam  extensional  modes  compared  to  bending 
modes)  are  discussed. 

A  goal  of  the  SOA  method  development  presented  here  is  to  improve  the 
efficiency  of  design  cycles  and  contribute  to  the  development  of  further  design 
optimization  tools,  especially  in  cases  where  structural  acoustic-enclosures  are  of 
interest  such  as  the  design  of  passenger  compartments  of  cars,  airplane  cabins, 
vibration  sensitive  payloads  of  space  vehicles  and  other  structures  where  noise 
and  vibration  levels  are  of  interest. 
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Chapter  1:  Order  Reduction  of  the  Coupled  Structural- Acoustic 

Problem 

1.1  Coupled  structural-acoustic  problem 

The  Equations  of  motion  for  an  acoustic  enclosure  coupled  to  a  flexible 
structure  is  written  as 

[MM+[K.lr}-pMAq,}=o  (» 

The  Equations  for  a  flexible  structure  with  such  an  enclosure,  including 
viscous  damping  and  an  applied  external  force,  is  given  by 

[M,fe,}  +  [C,te,}+fc]k,}-KrW=M  (2) 

where  the  matrices  [. Ma ]  and  [Ka]  are  the  acoustic  mass  and  stiffness  matrices 
respectively,  and  [a/j],[C5],  and  [*,],  are  the  structural  mass,  damping  and 
stiffness  matrices. 

The  matrix  [7?as]  is  the  structural-acoustic  coupling  matrix.  The  vectors 
{^  (/)},{/>(/)},  and  {F(0}  are  the  vectors  of  generalized  structural  displacement, 
pressures,  and  structural  excitation  forces.  p0  is  the  density  of  the  acoustic  fluid. 

Combining  equations  (1)  and  (2)  in  matrix  form,  the  coupled  structural- 


acoustic  equations  are  obtained  as 
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1.2  Full-order  structural-acoustic  model 


In  equation  (3),  the  coupled  structural-acoustic  mass  matrix  is 


(4) 


the  coupled  structural-acoustic  viscous-damping  matrix  is 


O' 

0 


(5) 


and  the  coupled  structural-acoustic  stiffness  matrix  is 


(6) 


The  Fourier-transformed  equations  for  a  structural  acoustic  system  in  the 
frequency  domain  can  then  be  written  as 

[-  co\M]  +  MC]  +  [K)\u(jco)}=  {F(jco)}  (7) 


where  [M\,  [C],  and  [Z]  are  the  matrices  as  shown  in  (4),  (5)  and  (6),  and 


{u{ja)}= 


fff.c/®)  1 

1  p(J<*>)  j 
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co  is  the  frequency  in  rad/s. 

Any  linear  responses  obtained  from  the  solution  of  equation  (7)  is 

y(jco)  =  {c}  *  {u(jco)}  =  {n(jco)}  *  {F(jco)}  (8) 

from  which  scalar  responses  such  as  stress,  displacement,  and  pressure  can  be 
calculated. 

The  full-order  adjoint  formulation  corresponding  to  equation  (7)  is 

(-  co2[M]  +  MC]  +  [K])*  {tj(j'co)}  =  {c}  (9) 

In  both  equations  (8)  and  (9),  the  superscript  *  denotes  the  complex 
conjugate  transpose  of  a  complex  matrix.  Note  that  unlike  in  pure  structural 
problems  where  the  mass  and  stiffness  matrices  [M\  and  [X]  are  symmetric,  in  the 
coupled  structural-acoustic  formulation  these  matrices  are  non-symmetric.  Also, 
in  the  adjoint  formulation  the  vector  {c}  is  frequency-independent. 

1.3  Reduced-order  structural-acoustic  model 


The  term  “full-order”  used  in  this  work  describes  the  complete  system  of 
equations  corresponding  to  a  finite-element  model  with  a  fine  mesh  and  a  large 
number  of  degrees  of  freedom.  A  “reduced-order”  system  is  then  a  system  of 


7 

equations  with  fewer  equations  used  to  approximate  the  “full-order”  system.  In 
the  Modal  Displacement  (MD)  method,  a  reduced-order  system  is  obtained  by 
multiplying  a  subset  of  modes  (shape  functions)  in  a  reduced-basis  matrix  [O] 
with  the  vector  of  generalized  displacements  {#(/©)}  as  written  below  in  equation 
(10). 

where  N  is  the  number  of  reduced-basis  vectors  used. 

Similarly  we  can  reduce  the  full-order  equations  given  in  (7).  If  we  write 
this  equation  as 

[AUco)]{u(j(o)}=  (fO©)}  (11) 

then  [A]  is  reduced  by  pre  and  post  multiplication  of  [<E>]*  and  [<f>]  such  that 

M*]=[«>]*M1[®]  (12) 

and  the  generalized  input  force  vector  {F}  is  also  reduced  such  that 
{Fo0©)}  =  [O]*{F(y®)}  (13) 

The  complete  reduced-order  direct  problem  is  then  written  as 


(14) 
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A  corresponding  reduced-order  adjoint  frequency  response  problem  can 
be  found  by  using  another  reduced  basis  matrix  [vp]  to  approximate  the  full-order 
adjoint  following  the  same  steps  used  for  the  direct  problem.  Starting  by  reducing 
the  generalized  adjoint  displacements  {r(ja>)},  such  that 

{f}  (0}  =  [fyq },  fo'# }]  K/®)}  =  P*']  {r(ja))  (15) 

Simplifying  the  equation  for  full-order  adjoint  equation  in  (9)  as 

[a]*{t]U(o)}={c}  (16) 

then  [A]*  is  reduced  by  pre  and  post  multiplication  by  the  basis  matrices  [vp]*  and 
[T]  such  that 

[A,]  =  m*Mm  (17) 

The  reduced  {c}  vector  is  then  given  by 

=  (18) 

The  complete  reduced-order  adjoint  problem  can  then  be  written  as 


[A^]*{r{j(o)}  =  {c^} 


(19) 


1.4  Selection  of  reduced-order  basis  matrices 
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The  reduced-order  basis  matrices  [<£>]  and  ['P]  are  selected  while  taking 
into  account  such  variables  as  spatial  distribution  of  excitation  forces,  weights 
used  to  determine  particular  responses  (the  vector  {c}),  frequency  oscillation  as 
well  as  material  properties.  Utilization  of  natural  modes  of  structures,  “fictitious 
mass”  modes,  Ritz  vectors  and  other  sub-component  modes  have  been  studied 
over  the  years  in  the  field  of  structural-dynamics  (Refs.  2,  8,  18,  19,  21,  26,  30, 
31,  35,  37).  In  the  general  structural-acoustic  cases  presented  here,  mode  shape 
matrices  (eigenvectors)  can  be  complex. 

In  this  work,  the  notation  [3>]  is  used  to  specifically  denote  the  reduced- 
order  basis  matrix  created  using  the  Modal  Displacement  (MD)  method  while  the 
notation  [¥]  is  used  to  specifically  denote  the  reduced-order  basis  matrix  created 
using  the  Ritz  Vector  (RV)  approach. 

It  is  also  noteworthy  that  [<£]  and  [¥]  can  be  interchanged  between  either 
the  direct  or  adjoint  formulation  as  long  as  it  is  used  consistently.  In  other  words, 
the  direct  problem  can  be  reduced  by  [<f>]  and  the  adjoint  problem  reduced  by  [VP], 
or  the  direct  problem  can  be  reduced  by  [T]  and  the  adjoint  problem  reduced  by 
[<t>].  The  advantages  and  shortcomings  of  either  approach  when  it  comes  to 
creating  first-order  approximations  to  be  used  for  the  Second  Order 
Approximation  (SOA)  will  be  presented  and  discussed  in  following  chapters. 
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Chapter  2:  A  Second-Order  Approximation  using  Combined 
Direct  and  Adjoint  Solutions 

2.1  Second-order  approximation  of  scalar  responses 

A  linear  problem  represented  as  a  system  of  equations  in  matrix  form 

MM-W  <2°) 

can  be  solved  for  {*}  directly: 

{*}  =  W,{6}  '  (21) 

This  system  also  has  the  complementary  adjoint  problem 

KrM=M  (22) 

A  full-order  scalar  response  y  can  be  found  using  either  the  direct  or  the  adjoint 
solution,  such  that 

y  =  {c}*  {*}  =  {c}  *  [A]'1  {b}  =  {17}  *  {b}  (23) 

Approximate  solutions  of  the  direct  (20)  and  adjoint  (22)  problems  that 
can  be  obtained  by  solving  equations  (10)  and  (14)  or  (15)  and  (19),  and  the 
reduced-order  scalar  response  yean  be  written  as 


y  =  {<?}*  {*}+ }*  M- {^7  }*  [^K?} 


(24) 
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The  ~  (tilde)  symbol  here  is  used  to  indicate  that  a  quantity  is  an  approximation 
of  the  variable  named.  Variables  that  do  not  have  the  ~  (tilde)  symbol  are  full- 
order.  The  errors  due  to  these  approximations  can  be  expressed  as 

{*}  =  {*}“{*} 

W}  =  M-{>7}  (25) 

5y  =  y-y 

or 

{*}={*}“{*} 

M={ y}-{sn}  (26) 

y=y-Sy 

Substituting  these  values  into  equation  (24)  shows  that  the  error  is  second  order 

Sy  =  -Sr]  *  A  Sx  (27) 

For  a  discussion  on  the  second-order  nature  of  this  approximation,  please 
consult  Appendix  D. 

2.2  Second  order  approximation  of  sensitivities 

Design  sensitivities  are  obtained  by  taking  the  derivative  of  the  scalar 
response  of  interest  with  respect  to  the  design  variable  of  interest.  For  the 
reduced-order  scalar  response  as  written  in  equation  (24)  then,  differentiating  with 
respect  to  the  design  variable  DV  gives 


12 


( W  *  P)  +  m  *  {»)  -  ®  *  M  {*))  (28) 

oDV  oDV 


Here,  the  sensitivities  of  the  approximate  direct  and  approximate  adjoint 
solutions,  8{x}/8DV  and  8{rf}/dDV  are  used  along  with  the  sensitivities  of  the 
full-order  vectors  and  matrices 8{b}/8DV ,8{c}/8DV and  d[A\l 8DV  to  write  the 
second  order  approximation. 

Expanding  equation  (28), 


dy 


8DV  8DV 
8{7j}* 


8{c}  *  „  ,  .  .  9{x}  8{tj}  *  * 


8DV 


8DV 

8[A] 


8DV 


m 

8DV 


8DV 


8{x}_ 

8DV 


(29) 


and  collecting  terms. 


=  ({c}  *  *  [A])^-  +  ({&}  -  [A]  {x }) 

8DV  V  8DV  8DV  V  ’ 


-{nr 


^}+^}+{ffym 


8DV 


8DV 


8DV 


(30) 


The  term  ({f»}-[^]{x})  is  the  error  vector  of  the  full  order  equation  (20) 
when  the  approximate  direct  solution  is  used  to  replace  the  exact  solution. 
Likewise  the  term  ({c}  *  -{77 }  *  [yf])  is  the  error  vector  of  the  full  order  equation 
(22)  when  the  approximate  adjoint  solution  is  used  instead  of  the  exact  solution.  If 
the  approximate  direct  and  approximate  adjoint  solutions  are  accurate,  that  is  the 
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quantities  ({6}-|>l]{x})  and  ({c}*-{^}*[^])  are  small,  then  the  equation  for 
the  sensitivity  of  the  second  order  approximation  can  be  simplified  as 


-{?>*- 
dDV  dDV 


d[A]  5{c}* 


dDV 


dDV 


(31) 


Furthermore,  recognizing  that  in  many  practical  design  optimization 
problems  the  {b}  and  {c}  vectors  remain  independent  of  the  design  variables,  the 
equation  in  (31)  can  be  simplified  by  realizing  that  the  derivatives 


d{c}* 

dDV 


(32) 


and 


M.o 

8DV 


Thus,  equation  (31)  simplifies  to 


& 

dDV 


*  *  M 


dDV 


{*} 


(34) 


In  certain  situations  where  {c}  does  depend  on  design  variables,  for 
instance  in  the  case  of  stress  constraints  where  {c}  depends  on  the  element  cross- 
sectional  area  parameter,  a  simplified  derivative  of  the  2nd  order  approximate 
behavior  output  is 
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dDV 


i-{7} 


dDV  dDV 


M! 

dDV 


-&}' 


d[A] 

dDV 


{*}  (35) 


2.3  Second  order  approximations  for  the  structural-acoustic  problem 


Applications  of  Order-reduction  techniques  for  structural-acoustic 
problems  are  widely  in  use  today  (Refs.  26,  29,  35).  In  most  applications,  natural 
modes  of  the  structure  are  used  in  the  order  reduction  of  the  structural  model. 
These  are  then  coupled  with  a  reduced-order  acoustic  model  based  on 
superposition  of  some  acoustic  modes  in  a  variant  of  the  substructure  modal 
synthesis  method  of  Structural  Dynamics  (Ref.  8).  Such  approaches  have  often 
shown  accurate  structural-displacement  responses  but  failed  to  provide  accurate 
acoustic  pressure  responses  (Ref  24).  This  is  because  in  structural-acoustic 
systems,  natural  modes  of  the  acoustic  system  are  normally  in  a  much  higher 
frequency  range  than  those  of  the  structural  part  of  the  model  (Refs.  10,  40).  For 
example,  it  is  desired  to  create  a  basis  matrix  from  the  first  40  modes  of  a  system 
to  analyze  a  system’s  behavior  over  a  range  of  frequencies.  In  the  case  of  a 
coupled  structural-acoustic  system,  let  the  40th  mode  be  at  100  Hz.  However,  the 
first  acoustic  mode  is  at  120  Hz.  Clearly  then  the  basis  matrix  will  only  contain 
structural  modes  and  no  acoustic  modes,  resulting  in  a  very  poor  approximation  of 
any  acoustic  results.  Even  if  more  modes  are  used,  the  likely  end  result  is  that  still 
a  majority  of  modes  selected  are  structural  modes  and  only  a  few  acoustic  modes 
are  included,  leading  to  similarly  poor  approximations  as  before.  Furthermore, 
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using  more  modes  also  means  that  the  approximation  is  less  efficient  as  more 
and  more  equations  are  added  to  the  system  of  equations. 

A  method  proposed  here  to  get  around  this  problem  is  by  creating  a 
pseudo  model  for  which  the  structural  and  acoustic  modes  are  closer  together  to 
allow  more  structural  and  acoustic  modes  to  be  used  when  creating  basis- 
matrices.  This  is  achieved  by  altering  material  properties,  such  as  fluid  density,  to 
shift  the  frequency  at  which  natural  modes  occur.  Naturally,  this  pseudo  model  is 
only  used  for  creating  basis  matrices,  and  not  for  the  analysis  itself.  Increasing 
density  of  an  acoustic  fluid  means  that  the  fluids  modal  frequencies  move  lower 
in  the  frequency  band  towards  where  the  structural  modes  reside*.  Now  the  basis 
matrix  created  will  include  a  significant  number  of  both  structural  and  acoustic 
mode  shapes.  A  comparison  for  accuracy  between  approximations  using  this 
method  and  one  where  this  aspect  is  ignored  is  included  later  in  this  thesis. 

As  a  final  remark,  for  the  structural-acoustic  coupled  problem  when 
modes  to  create  the  basis  matrix  [<£>]  are  used,  these  are  then  in  the  present  work 
the  coupled  modes  of  the  structural-acoustic  system.  When  Ritz  vectors  are  used 
to  create  the  basis  matrix  [¥]  they  are  made  of  responses  of  the  coupled  full-order 
structural-acoustic  system  to  various  inputs. 

♦Note:  By  “structural  modes”  and  “acoustic  modes”  in  the  case  of  a 
coupled  system,  we  mean  coupled  modes  in  which  either  the  structural 
contributions  or  acoustic  contributions  are  dominant,  respectively. 
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Chapter  3:  The  Design-Oriented  Structural- Acoustic  2D  Finite 

Element  Model 

3.1  Finite-element  formulation 

The  design-oriented  structural-acoustics  finite  element  model  used  here 
for  2D  problems  is  based  on  beam  finite  elements  for  the  structural  boundaries  of 
an  acoustic  enclosure,  coupled  with  triangular,  linear-pressure  elements  for  the 
acoustic  fluid.  A  subroutine  was  written  in  Fortran  90  to  create  the  necessary 
matrices,  both  the  coupled  structural-acoustic  stiffness  [AT]  and  mass  matrices  [M\, 
and  also  the  sensitivity  matrices  with  respect  to  the  design  variables  of  interest. 
The  design  variables  of  interest  are  the  cross-sectional  areas  of  the  beams  (A)  and 
the  moment  of  inertia  (I).  Hence  the  sensitivity  matrices  created  are  [dK  /  dA], 
\dK  I dl\,  \dM  I dA\,  and  [dM/dl].  Detailed  derivations  for  the  acoustic 
equations  are  presented  in  Appendix  A,  while  those  for  the  acoustic  component  of 
the  mass  and  stiffness  matrices  are  presented  in  Appendix  B.  Appendix  C  outlines 
the  structural-acoustic  interaction  terms,  as  well  as  the  combined  structural- 
acoustic  formulation. 

After  these  matrices  were  created,  they  were  imported  into  MATLAB, 
which  facilitates  easier  matrix  manipulations,  for  calculation  of  dynamic 
deformations,  internal  stresses  in  structural  elements,  pressures  at  node  on  the 
acoustic  system.  Sensitivities  with  respect  to  the  above  mentioned  design 
variables  are  also  calculated.  Direct  and  adjoint  problem  solutions  are  available  as 
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well,  both  in  full-order  and  approximate  order  using  the  various  order  reduction 
methods  described. 

3.2  Description  of  the  finite  element  model 

The  test  case  is  a  2-dimensional  acoustic  enclosure  surrounded  by  flexible 
walls  on  6  sides  as  shown  in  Figure  1.  The  model  is  inspired  by  previous  work  in 
structural-acoustic  finite  element  analysis  (Ref.  33). 

For  each  analysis,  a  reference  structure  and  a  modified  structure  are 
created.  The  reference  structure  is  used  to  obtain  the  basis  matrices  needed  to 
create  a  reduced-order  model.  The  modified  structure  is  one  where  the  element 
labeled  “Variable  Element”  is  modified  due  to  some  changes  made  to  its  design 
variables. 

The  modified  structure  is  then  solved  in  full  order  to  establish  a 
benchmark  solution.  This  benchmark  is  then  used  to  evaluate  accuracy  of  the 
various  approximations  obtained  from  order-reduced  methods.  In  design 
optimization  problems,  the  basis  matrices  only  have  to  be  obtained  once  from  the 
reference  structure,  and  then  can  be  used  to  approximate  solutions  of  many 
modified  structures  for  successive  changes  in  design  variables,  resulting  in 
significant  computational  time  savings. 


3.2.1  Geometry  of  the  finite-element  model 
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The  model  shown  in  Figure  1  has  143  acoustic  nodes,  44  structural  nodes, 
and  a  total  of  275  degrees  of  freedom.  An  additional  model,  identical  in  geometry 
but  using  a  much  finer  finite  element  mesh  was  also  created  and  is  shown  in 
Figure  2.  This  model  contains  1147  acoustical  nodes,  132  structural  nodes  and  a 
total  of  1 543  degrees  of  freedom,  and  will  be  referred  to  as  the  large-scale  model 
throughout  this  work.  This  large-scale  model  is  used  to  demonstrate  the 
computational  time  saving  achieved  by  reducing  a  1543  x  1543  system  to  a  much 
smaller  and  more  manageable  system. 

As  an  illustration  for  the  time  saving  that  can  be  realized,  a  test  case  of  a 
large-scale  full  order  system  (1543  x  1543  system)  requires  about  1  hour  and  40 
minutes  to  solve  on  a  computer  with  a  2.66  GHz  processor  with  512  MB  of  RAM, 
for  one  iteration  over  200  steps  in  the  frequency  domain.  In  order  to  solve  this 
system  using  the  second-order  method  utilizing  Ritz  vectors  and  Mode 
displacements  as  the  first  order  approximations,  an  initial  run  is  necessary  to 
obtain  the  basis  matrices.  This  initial  run  takes  approximately  2  hours  and  20 
minutes.  However,  after  this  initial  investment,  the  system  will  be  reduced  and 
each  additional  run  only  requires  2.5  minutes  (in  the  case  of  a  system  reduced  to  a 
40  x  40  system). 

Both  small-scale  and  large-scale  models  represent  a  2  dimensional,  six- 
sided  acoustic  enclosure,  where  the  acoustic  fluid  is  completely  enclosed  by  a 
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flexible  beam  wall,  and  is  supported  by  pivots  in  3  comers.  The  enclosure  is 
0.508  m  wide  on  the  top  and  bottom,  and  0.8763  m  high.  If  the  lower  left  comer 
of  the  enclosure  is  labeled  vertex  1  and  is  located  at  the  origin,  then  going 
counter-clockwise,  the  coordinates  of  the  vertices  are  as  outlined  in  Table  1 . 


Table  1:  Vertex  coordinates  of  acoustic  enclosure 


Vertex 

x-coordinate  y- 

coordinate 

(m) 

(m) 

1 

0 

0 

2 

0.508 

0 

3 

0.8453 

0.5842 

4 

0.508 

1.1684 

5 

0 

1.1684 

6 

-0.506 

0.8763 

The  beams  connecting  the  vertices  and  surrounding  the  enclosure  are 
constrained  in  the  x-y  direction  at  vertices  1,  2  and  4,  but  otherwise  are  allowed  to 
rotate  freely.  All  other  vertices  are  allowed  to  move  freely.  Also,  similar  to  the 
model  in  Ref.  33,  two  point  masses  of  0.03  kg  were  attached  to  the  bottom  beam 
panel  at  jc  =  0.2032  m  and  0.3048  m.  The  excitation  force  is  a  point  unit  load 
located  at  x  =  0.3048  m. 

3.2.3  Material  properties  of  the  reference  model 

The  beams  surrounding  the  enclosure  are  made  of  aluminum,  with 
modulus  of  elasticity  (E)  of  72  GPa  and  density  of  2600  kg/m3.  The  beams  have 
rectangular  cross  sections  with  0.2032  m  at  the  base,  and  height  of  0.003175  m. 
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Thus  they  have  a  cross-sectional  area  of  6.45  x  10‘4  m2  and  a  moment  of  inertia 
of  5.4197  x  10'10  m4.  The  fluid  in  the  enclosure  is  air,  with  density  of  1.225  kg/m3, 
and  speed  of  sound  of  342  m/s. 

3.2.4  Material  properties  of  the  modified  models 

The  modified  models  all  maintained  the  same  geometrical  properties  as 
the  reference  model.  However,  in  the  element  marked  as  “Variable  Element”  in 
Figure  1,  the  design  variables  of  interest  were  changed.  Both  the  moment  of 
inertia  and  beam  cross-sectional  area  were  increased  by  50%,  corresponding  to  I = 
8.1296  x  10'10  m4  and  A  =  9.675  x  10-4  m2.  This  model  was  created  to  represent  a 
design  optimization  problem  where  the  design  variables  are  changed  as  the  model 
is  driven  to  the  optimum. 

Several  models  were  also  created  with  changes  in  the  design  variables  I 
and  A  varying  between  10%  to  100%,  and  are  used  to  illustrate  how  the 
approximations  behave  as  the  modified  model  moves  further  away  from  the 
reference  model. 

Finally,  two  different  values  of  the  structural  damping  factor  (g)  are  used, 
0.06  and  0.15.  While  the  value  of  the  damping  factor  itself  does  not  affect  the 
accuracy  of  the  approximate  solution  in  a  fundamental  way,  a  higher  damping  has 
the  effect  of  smoothing  out  sharp  peaks  in  the  response  curves,  which  in  turn 
reduces  the  difference  in  error  calculations  caused  by  small  shifts  of  very  narrow 
frequency  response  peaks. 


21 

Another  model  that  was  created  is  what  will  be  referred  to  as  the  pseudo 
model.  This  model  is  created  because  the  difficulties  of  modal  approximation 
methods  in  finding  proper  reduced  bases  for  reducing  the  order  of  structural- 
acoustic  models  due  to  the  large  differences  in  modal  frequencies  between  the 
structural  and  the  acoustic  modes  (Refs  10,  24,  40).  This  model  is  used  only  to 
obtain  the  modal  basis  matrices  for  the  modal  approximation  methods,  as  it  does 
not  have  any  physical  meaning,  other  than  that  it  is  a  model  where  the  structural 
and  acoustic  modes  are  closer  together  in  the  frequency  range,  while  the 
corresponding  mode  shapes  remain  the  same.  To  achieve  this,  the  fluid  density  in 
this  model  was  increased  100  times,  and  the  cross-sectional  areas  of  the  beams 
were  reduced  to  0.01  times  of  the  original.  (Please  see  note  on  page  15). 

3.3  Validation  of  the  finite  element  model 

To  validate  the  model  used  in  this  work  with  dedicated  codes,  two 
comparison  tests  were  made  using  a  similar  model  analyzed  by  the  commercial 
FEM  software  ANSYS  (Ref.  1).  The  first  test  is  a  comparison  of  natural 
frequencies  between  the  two  models.  Table  2,  which  records  the  first  12  natural 
modes  of  both  systems,  shows  that  the  natural  modes  of  the  two  systems  are 
identical,  as  expected. 


Table  2:  Comparison  of  natural  frequencies. 
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Mode 

ANSYS  FEM  Code 

1 

(Hz) 

6.7 

(Hz) 

6.7 

2 

18.0 

18.1 

3 

21.4 

21.4 

4 

28.2 

28.3 

5 

31.9 

32.0 

6 

40.4 

40.3 

7 

50.2 

50.4 

8 

68.8 

68.0 

9 

73.5 

73.4 

10 

85.5 

85.6 

11 

91.2 

91.2 

12 

107.0 

107.3 

The  second  test  is  a  comparison  of  sound  pressure  output  in  decibels  at  the 
node  marked  in  Figure  1.  The  result  again  shows  very  good  correlation  between 
the  two  models  and  is  presented  in  Figure  2.  The  sound  pressure  output  calculated 
from  the  ANSYS  model  are  virtually  identical  to  that  calculated  by  the  FEM 
program  written  for  this  work  over  the  frequency  range  of  25  Hz  to  250  Hz.  Based 
on  these  two  tests,  it  is  shown  that  the  FEM  model  created  for  this  work  is  indeed 


accurate. 


Figure  1:  2D  Structural  Acoustic  System 
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Figure  2:  Large-scale  2D  Structural  Acoustic  System 


Chapter  4:  Numerical  Results 


26 


4.1  Introduction  to  the  test  cases 

The  Second  Order  Approximation  results  were  obtained  using  several  different 
approaches. 

In  the  first  approach,  the  Ritz  vector  is  used  to  reduce  the  order  of  the 
direct  problem  written  in  Eq.  7  corresponds  to  the  force  input  on  the  right  hand 
side  of  the  equation.  The  reduced  basis  [¥</]  is  also  frequency-dependent  and 
consists  of  single  frequency-dependent  vectors  each  a  solution  of  the  full  order 
reference  problem.  The  subscript  d  here  is  used  to  denote  that  this  basis  vector 
corresponds  to  the  direct  problem.  When  Ritz  vectors  were  used  to  reduce  the 
order  of  the  direct  problem,  then  the  adjoint  problem  was  reduced  by  a  subset  of 
the  left  natural  modes  [<!></]. 

(36) 

In  the  second  approach,  a  Ritz  vector  approximation  is  used  for  the  adjoint 
problem  corresponds  to  a  response  y(jco)  as  given  in  Eq.  8.  For  this,  the  reduced 
basis  ['Fa]  used  is  the  solution  of  Eq.  9.  This  basis  [TJ  is  frequency  dependent, 
and  varies  with  frequency.  The  subscript  a  here  is  used  to  denote  that  this  basis 
vector  corresponds  to  the  adjoint  problem.  The  reduced-order  Ritz  solution  is  of 
order  1  x  1  based  on  a  single  adjoint  solution  at  each  frequency  of  the  reference 
system.  Since  Ritz  vectors  were  used  to  reduce  the  order  of  the  adjoint  problem. 
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then  a  low-frequency  subset  of  the  right  natural  modes  [<Da]  was  used  to  reduce 
the  direct  problem. 

(Ai[M)  +  [KM)  =  {  0}  (37) 

Since  the  structural-acoustic  mass  and  stiffness  matrices  in  the  formulation 
used  here  are  non-symmetric,  the  mode  shapes  are  generally  complex.  Symmetric 
formulations  for  mass  and  stiffness  matrices  (Refs  9,  29)  can  be  used,  however 
none  of  the  resulting  matrices  are  positive  definite,  hence  the  modes  will 
generally  remain  complex. 

In  another  approach,  the  adjoint  and  direct  approximations  were  using  the 
adjoint  and  direct  solutions  of  the  original  reference  model.  This  saves 
computational  time,  as  there  is  no  need  to  calculate  the  first  order  approximations 
using  reduced  order  models.  A  variation  of  this  approach  is  to  use  one  of  the 
original  reference  solutions  and  combine  it  with  a  first  order  approximation  (such 
as  the  modal  approximation)  to  improve  results. 

Finally,  first  order  Taylor  series  approximations  were  used  to  create  the 
approximate  first  order  results.  In  this  approach,  the  approximate  first  order  basis 
for  the  second  order  approximation  takes  the  form 
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In  this  case,  the  sensitivity  derivative  -  is  the  sensitivity  of  the 

dDV,  0 

original  reference  model.  In  the  following  sections,  second-order  approximation 
results  obtained  using  various  first-order  approximations  will  be  presented  and 
evaluated.  In  all  the  studies  presented  here,  the  following  steps  were  taken.  First, 
analysis  and  sensitivity  results  were  generated  for  the  reference  model,  and 
reduced  basis  vectors/matrices  obtained.  Then  the  model  was  modified  by 
changing  properties  of  the  variable  element.  A  full-order  detailed  analysis  of  the 
modified  model  was  carried  out  and  used  to  check  accuracy  of  the  approximate 
analysis  applied  to  this  modified  model.  Finally,  the  modified  model  was 
analyzed  using  a  variety  of  reduced-order  techniques  and  accuracy  of  the  reduced 
order  calculations  assessed.  In  an  optimization  environment  detailed  analyses  are 
carried  out  only  a  small  number  of  times  while  most  of  the  optimization  search 
effort  is  based  on  approximate  analyses. 

As  can  be  seen  a  large  number  of  combinations  of  first-order  methods 
were  used  to  obtain  the  second-order  approximation.  Some  combinations  are  not 
presented  here  as  they  did  not  add  additional  insight,  or  were  simply  inaccurate. 
Those  that  are  presented  are  summarized  in  Table  3,  and  are  discussed  further  in 


the  sections  that  follow. 


Table  3:  Summary  of  approximation  methods 
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Case 

First-order  Direct 

First-order  Adjoint 

Scale  of  Model 

i 

Mode  Displacement 

Ritz  Vectors 

Small 

2 

Ritz  Vectors 

Mode  Displacement 

Small 

3 

Original  Solution 

Original  Solution 

Small 

4 

Mode  Displacement 

Original  Solution 

Small 

5 

First-order  Taylor 

First-order  Taylor 

Small 

6 

Mode  Displacement 

Ritz  Vectors 

Large 

7 

Ritz  Vectors 

Mode  Displacement 

Large 

4.2  Dynamic  results  of  the  second  order  approximation,  small  scale  model 
4.2.1  Adjoint  based  on  Ritz  Vectors,  Direct  Based  on  Modes 

In  Figure  4  through  Figure  9,  the  reduced  order  direct  solution  is  based  on 
a  single  Ritz  vector  (1  x  1),  and  the  adjoint  solution  is  based  on  a  reduced  basis 
matrix  of  41  modes  (41  x  41).  These  numbers  of  modes  chosen  were  such  that  it 
covers  at  least  2  times  the  frequency  range  of  interest  in  the  original 
configuration. 

The  structure  is  excited  harmonically  at  the  bottom  wall,  and  errors  of  the 
reduced  approximation  presented  versus  the  frequency  of  excitation,  between  20 
Hz  and  260  Hz. 

As  can  be  seen,  the  Ritz  approximation  has  the  least  accuracy,  relying  on 
only  1  Ritz  vector  to  capture  the  behavior  of  the  system,  since  there  is  only  1 


excitation  force. 
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The  modal  approximation  shows  better  accuracy,  especially  in  the  lower 
ranges  of  the  frequency  band,  but  starts  to  diverge  above  the  200  Hz  range.  This 
can  be  attributed  to  the  lack  of  higher  frequency  modes  in  the  basis  matrix  due  to 
only  using  the  first  41  modes.  Improvements  can  be  achieved  by  adding  more 
modes  to  the  basis  matrix,  but  at  an  increase  of  computational  cost. 

By  combining  the  1  x  1  Ritz  approximation  with  the  40  x  40  modal 
approximation,  the  second-order  method  results  show  marked  improvement  of 
accuracy.  It  is  important  to  note  that  while  the  modal  approximation  solution 
plotted  was  approximated  using  41  modes,  the  second  order  approximation  only 
used  40  modes,  one  less  mode  in  order  to  give  a  fair  comparison  since  the  second 
order  method  “gains”  one  extra  approximation  vector  from  the  Ritz 
approximation.  This  is  done  to  show  that  the  second-order  approximation  does  not 
gain  extra  accuracy  merely  from  using  more  basis  vectors. 

In  Figure  4  errors  in  the  approximate  bending  moment  (stress)  are  shown 
for  the  element  shown  in  Figure  1.  Figure  5  shows  the  error  of  the  different 
approximation  methods  compared  to  the  full-order  solution  when  computing 
pressure  at  the  node  indicated  in  Figure  1.  Again  the  reduced-order  direct  solution 
is  based  on  a  single  Ritz  vector  (1  x  1),  and  the  adjoint  solution  is  based  on  a 
reduced  basis  matrix  of  41  modes  (41  x  41).  The  model  is  excited  with  the  same 
load.  The  Ritz  vector  approximation  was  generally  not  as  accurate  as  the  modal 
approximation,  while  the  second-order  approach  showed  improvement  in 
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accuracy,  particularly  in  the  upper  frequency  range  where  the  modal 
approximation  begins  to  diverge. 

Accuracy  of  predicting  reduced  order  design  sensitivities  are  examined 
/  next.  Figure  6  shows  errors  in  the  different  methods  of  approximation  relative  to 
the  full-order  solution  when  predicting  sensitivities  of  bending  moment  with 
respect  to  changes  in  the  beam  moment  of  inertia  ( dM/dl ).  Here  is  where  some 
significant  errors  in  the  first-order  methods  are  observable.  The  error  in  the 

t 

approximate  direct  solution  occasionally  exceeds  60%,  while  the  error  for  the 
approximate  adjoint  solution  exceeds  30%  at  some  point.  By  contrast,  the  second 
order  approximation  never  exceeds  5%  in  the  frequency  range  shown. 

In  Figure  7,  errors  in  the  different  methods  of  approximation  when 
predicting  sensitivities  of  bending  moment  with  respect  to  changes  in  the  beam 
cross-sectional  area  are  presented.  It  becomes  evident  here  that  errors  with  the 
direct  approximation  are  uniformly  quite  large.  With  only  a  few  exceptions,  the 
error  of  this  method  routinely  exceeds  40%.  In  two  areas,  between  30  Hz  and  40 
Hz,  and  between  170  Hz  and  175  Hz,  the  error  even  exceeds  500%.  The  adjoint 
problem  approximation  behaved  better  with  errors  never  exceeding  25%.  The 
second-order  method  however  only  performs  marginally  better  here.  Since  it  uses 
information  from  two  first-order  approximation  methods,  if  one  of  them  is 
significantly  inaccurate,  the  second-order  method  will  typically  only  perform 
marginally  better  or  equivalent  to  the  more  accurate  of  the  two.  This  pattern  is 
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repeated  in  all  results  where  one  of  the  first-order  approximations  has  errors  far 
exceeding  100%. 

Sensitivities  of  pressure  with  respect  to  design  variables  were  also 
examined.  In  Figure  8,  errors  in  the  different  methods  of  approximation  when 
predicting  sensitivities  of  pressure  with  respect  to  changes  in  the  moment  of 
inertia  ( dP/dl)  is  shown.  As  was  observed  in  the  dMJdA  results,  the  direct 
approximation  solution  showed  errors  larger  than  150%  throughout  the  frequency 
range.  The  adjoint  approximation  solution  also  shows  errors  of  up  to  150%  in 
some  areas.  As  noted  previously,  if  one  of  the  first-order  approximations  shows 
large  errors,  then  the  second-order  approximation  method  will  only  perform  as 
good  as  the  better  first-order  approximation.  Still,  throughout  the  frequency  range, 
the  maximum  error  of  the  second-order  approximation  never  exceeds  100%, 
while  both  first-order  approximations  had  a  maximum  error  of  150%  in  the  same 
frequency  range. 

In  Figure  9  we  examine  errors  of  the  different  approximation  methods 
when  predicting  the  sensitivity  of  the  pressure  with  respect  to  the  beam  cross- 
sectional  area  (dP/dA).  Again,  despite  the  large  errors  in  the  predictions  of  the 
direct  approximation,  the  second-order  solution  generally  is  at  least  as  good  as  the 
adjoint  approximation.  The  exception  being  in  the  regions  above  200  Hz  where 
the  second-order  solution  generally  had  less  than  half  of  the  error  of  the  adjoint 
approximation,  and  significantly  more  accurate  than  the  direct  approximation 
using  Ritz  vectors. 
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4.2.2  Adjoint  based  on  modes,  direct  based  on  Ritz  vectors 

In  this  approach,  a  Ritz  vector  approximation  was  used  for  the  adjoint 
problem,  while  the  modal  approximation  was  used  for  the  direct  problem.  Here 
the  Ritz  vector  corresponds  to  a  response  y(j<a).  That  is,  each  response  has  its  own 
{c}  vector  (Eq.  19)  and  it  is  used  as  a  right-hand  side  in  the  adjoint  equation  to 
produce  the  Ritz  vector.  In  Figure  10,  it  is  shown  how  this  approach  performs 
when  approximating  the  bending  response.  As  in  previous  cases,  the  second-order 
bending  response  approximation  has  consistently  smaller  errors  when  compared 
to  the  first-order  approximations. 

However,  a  problem  that  is  observed  when  using  this  method  is  that  the 
Ritz  vector  approximation  is  simply  unable  to  capture  the  pressure  and  pressure 
related  derivatives.  This  could  be  due  to  the  response  at  the  acoustic  element  is 
relatively  an  order  of  magnitude  smaller  when  compared  to  the  response  at  the 
beam  element.  This  can  be  observed  in  Figure  1 1.  In  this  approximation,  the  error 
due  to  the  adjoint  based  Ritz  vector  approximation  is  so  large  it  is  outside  of  the 
chart.  In  cases  where  one  first-order  approximation  is  .  highly  inaccurate,  the 
second-order  method  performs  only  as  well  as  the  better  first-order  method,  in 
this  case  the  modal  approximation. 


4.2.3  Adjoint  and  direct  based  on  original  solution 
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Here,  the  original  solution,  i.e.  solution  of  the  unmodified  reference 
model,  was  used  as  the  first-order  approximation  for  the  second-order  method. 
Also,  changes  in  the  design  variables  of  the  variable  element  were  increased  to 
100%,  and  the  structural  damping  factor  g  is  reduced  to  0.06.  Figure  12  shows  the 
pressure  response  for  the  entire  range  of  frequencies,  while  Figure  13  shows  the 
same  for  a  zoomed  area.  It  can  be  seen  here  that  the  second-order  response 
approaches  the  modified  full-order  solution  is  only  somewhat  closer  than  just  the 
original  full-order  solution  alone.  Examining  the  error  plot  in  Figure  14  confirms 
that  the  second-order  approximation  shows  some  improvement  over  the  original 
solution. 

Examining  the  sensitivities,  Figures  15  and  16  show  the  errors  of  the 
approximate  solutions  when  estimating  the  sensitivity  of  the  moment  and  pressure 
responses  with  respect  to  the  cross-section  Area  (A)  design  variable.  In  this  case, 
error  of  the  second-order  method  is  in  some  cases  smaller  by  almost  a  factor  of  4 
or  more  (25%  vs.  100%)  in  various  regions  of  the  frequency  domain. 

4.2.4  Adjoint  based  on  modes,  direct  based  on  original  solution 

Here  it  is  sought  to  improve  on  the  previous  method  by  using  a  more 
accurate  first-order  basis.  Having  performed  well  in  the  previous  approximation 
methods,  the  modal  basis  approximation  is  used  to  replace  one  of  the  first-order 
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original  solutions.  The  result  is  a  method  that  combines  the  adjoint  solution 
approximation  using  modal  basis  vectors  with  an  approximation  using  the  direct 
solution  of  the  original  (unmodified)  model.  These  two  approximations  are  then 
used  to  obtain  the  second-order  approximation.  From  a  computational  cost 
perspective  this  method  requires  more  resources,  as  the  modal  basis  matrix  has  to 
be  calculated  first.  However,  as  in  previous  methods,  once  this  basis  matrix  is 
calculated,  it  can  be  reused  for  future  iterations  for  all  changes  in  the  design 
variables. 

For  the  case  examined,  the  modified  beam  cross-sectional  area  and 
moment  of  inertia  were  increased  100%  from  the  reference  model,  and  a  damping 
factor  g=0.06  was  used.  Figure  17  shows  the  pressure  response  around  the  same 
region  as  shown  in  Figure  13.  Comparing  these  responses  shows  that  the  second- 
order  result  is  now  improved  over  the  results  previously  shown  in  Figure  13. 
Examining  the  error  of  the  second-order  pressure  response  as  shown  in  Figure  18 
confirms  that  over  the  entire  frequency  range  the  second-order  method  gives 
significant  improvement  over  the  first-order  methods.  Similarly  Figures  19  and  20 
show  improved  second-order  approximation  results  for  bending  sensitivities, 
while  Figures  21  and  22  show  error  of  the  second-order  approximation  results  for 
pressure  sensitivities.  In  all  these  sensitivity  results,  the  second-order 
approximation  in  general  shows  less  error  than  any  of  the  first-order  methods 
alone.  Also,  throughout  the  range,  error  of  the  second-order  method  is  more 
consistent.  For  example,  in  Figure  18,  around  200  Hz  -220  Hz,  the  first-order 
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approximations  exhibit  large  errors  to  over  80%,  while  the  second-order 
method  never  exceeds  20%  error  throughout  the  whole  frequency  band. 

.  As  noted  before,  the  second-order  method  results  improve  upon  the  first 
order  method  results,  and  when  results  from  the  first  order  methods  become  too 
inaccurate,  the  second  order  results  typically  also  become  less  accurate.  This  can 
be  observed  on  Figure  21,  particularly  in  the  range  between  200  Hz  and  250  Hz. 
As  the  error  in  the  first-order  methods  becomes  large,  the  second-order  error  too 
becomes  large.  However,  it  is  important  to  note  that  in  this  particular  case,  while 
the  second-order  error  is  relatively  large,  it  still  improves  on  the  first-order 
method  results. 

To  see  whether  the  approximation  results  can  be  better  captured  in  the 
case  of  models  with  smaller  changes  to  the  design  variables,  the  model  was  solved 
approximately  again  with  the  same  method,  except  the  design  variables  A  and  I 
were  only  changed  50%,  compared  to  100%  previously.  Also,  the  damping 
coefficient  was  increased  to  g  =  0.15  to  smooth  out  and  widen  the  resonance 
frequency  peaks.  The  desired  effect  of  this  is  to  reduce  errors  due  to  shifting  of 
the  peaks.  Figures  23  and  24  contrast  the  pressure  response  between  these  two 
models.  In  Figure  23,  the  pressure  response  is  shown  for  the  model  with  100% 
change  in  DV  and  0.06  damping,  and  Figure  24  shows  the  pressure  response  for 
the  model  with  50%  change  in  the  two  DVs  and  0.15  damping.  As  it  can  be  seen 
from  comparing  these  two  figures,  some  of  the  resonance  peaks  are  indeed 
flattened,  as  expected. 
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Figures  25  and  26  show  error  in  the  approximations  of  the  pressure  and 
bending  responses  respectively,  in  the  model  with  50%  change  in  the  two  DVs. 
As  it  can  be  seen  the  second-order  method  is  now  clearly  far  superior  than  the 
first-order  methods,  showing  errors  that  are  less  than  2%,  compared  to  up  to  12- 
20%  and  even  higher  for  the  first-order  methods. 

Examining  the  response  plots  also  confirms  that  the  second-order 
method’s  precision  is  better  than  that  of  each  of  the  first-order  methods.  Figure  27 
shows  a  segment  of  the  pressure  response,  and  Figure  28  shows  a  segment  of  the 
bending  response  curves.  In  both  figures  the  results  are  magnified  for  a  better 
comparison. 

Errors  in  approximating  the  sensitivity  of  the  design  variables  of  interest 
are  also  generated.  As  it  can  be  seen  in  Figures  29  and  30  for  the  bending  moment 
sensitivity  with  respect  to  the  design  variables  /  and  A  respectively,  the  second- 
order  approximation  produces  the  most  consistently  accurate  approximations.  The 
same  observations  can  be  made  for  Figures  3 1  and  32  where  the  errors  in  pressure 
sensitivity  approximations  are  presented.  Another  observation  that  can  be  made  is 
that  the  magnitude  of  the  errors  in  the  second  order  approximation  is  related  to  the 
errors  in  the  first  order  approximation.  For  example,  in  Figure  32  near  the  100  Hz 
range,  error  of  both  of  the  first-order  approximation  is  large,  so  the  error  in  the 
second-order  approximation  is  also  large,  even  though  it  is  still  smaller  than  each 
of  the  first  order  errors.  In  Figure  31  at  the  200  Hz  range,  error  in  the  first-order 
adjoint  approximation  is  large,  however  the  error  for  the  original  full-order 
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solution  is  not  large.  Error  in  the  second-order  approximation  is  higher  here 
than  at  any  other  point,  but  still  much  smaller  than  the  direct  first-order. 

4.2.5  Adjoint  and  direct  based  on  Taylor-series  approximation 

Another  attempt  at  using  additional  first-order  approximations  was  made 
using  Taylor  series  to  obtain  the  approximate  results.  The  modified  model  for  this 
case  had  only  one  design  variable  changed,  in  this  case  the  moment  of  inertia  (I) 
of  the  variable  element.  The  design  variable  was  also  changed  incrementally  from 
10%  to  50%  to  see  the  effect  of  the  larger  change  on  the  accuracy  of  the  second- 
order  approximation.  Figures  33,  34,  35  and  36  show  errors  in  pressure 
approximation  using  this  method,  as  the  change  in  design  variable  is  increased  by 
10%,  20%,  30%  and  50%  respectively.  As  it  can  be  seen,  the  second-order 
approximation  increasingly  becomes  less  accurate  as  the  change  is  increased. 
Approximation  errors  of  the  sensitivities  were  also  obtained,  and  in  Figures  37, 
38,  39  and  40  errors  of  sensitivity  approximation  of  the  pressure  with  respect  to 
the  moment  of  inertia  of  the  beams  (< dP/dl)  are  presented.  As  before,  the  figures 
correspond  respectively  to  10%,  20%,  30%  and  50%  changes  in  the  design 
variable.  Using  the  Taylor  expansion  method,  the  first-order  adjoint  and  direct 
approximations  for  the  sensitivities  are  identical.  Similar  to  previous  results, 
accuracy  of  the  second-order  approximation  while  very  good  at  small  changes  in 
design  variables,  become  steadily  worse  as  the  changes  are  increased.  When  the 
design  variable  was  changed  by  more  than  50%  the  second-order  solution  ceases 
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to  give  more  accurate  results  than  the  best  first-order  Taylor  expansion 
approximation.  The  reason  is  that  one  of  the  first-order  approximations  being 
highly  inaccurate. 

4.3  Dynamic  results  of  the  second-order  approximation,  large  scale  model 
4.3.1  Adjoint  based  on  modes,  direct  based  on  Ritz  vectors 

Approximation  results  for  the  large-scale  model  are  of  great  interest, 
because  they  test  the  viability  of  solving  large  set  of  equations  with  an 
approximation  based  on  a  much  smaller  sub-set  of  equations.  As  already 
mentioned,  the  large-scale  model  has  1543  degrees  of  freedom.  Also  for  this 
problem,  modes  for  the  modal  approximation  were  obtained  from  a  model  for 
which  the  density  of  the  fluid  was  increased  lOOx  and  the  cross-sectional  area  of 
the  beams  reduced  0.0 lx.  This  is  to  simulate  situations  where  in  the  process  of 
optimization  basis  matrices  for  the  approximations  are  obtained  from  a  model  that 
is  very  different  than  the  model  being  solved.  This  has  the  advantage  of  not 
having  to  recalculate  the  modes  every  time  the  model  changes. 

Figure  41  shows  error  of  the  different  approximation  methods  compared  to 
the  full-order  solution  when  computing  pressure  at  the  node  of  interest.  As  in  the 
case  for  the  small-scale  model,  the  reduced-order  direct  solution  is  based  on  a 
single  Ritz  vector  (1  x  1),  and  the  adjoint  solution  is  based  on  reduced  basis 
matrix  of  50  modes  (50  x  50).  The  difference,  of  course,  being  that  the  model 
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reduced  now  contains  1543  degrees  of  freedom,  corresponding  to  a  matrix  of 
size  1543  x  1543.  As  previously  in  the  small  model,  change  in  design  variables  I 
and  A  are  50%. 

As  it  can  be  seen  in  Figure  41,  error  of  the  second-order  approach  remains 
consistently  the  smallest  of  the  three  approximations.  Figure  42  shows  the  actual 
response  of  the  approximations  and  plots  the  full-order  result  along  side  the  3 
approximations  in  the  20  Hz  to  60  Hz  range.  As  can  be  seen,  the  second-order 
approximation  solution  traces  the  full-order  solution  very  closely. 

In  Figure  43,  the  bending  moment  approximation  errors  are  presented.  The 
second-order  approximation  is  far  more  accurate  than  each  of  the  first-order 
methods.  The  largest  error  recorded  for  the  second-order  method  in  the  entire 
frequency  range  did  not  exceed  20%  at  any  time. 

Figures  44  and  45  examine  the  bending  moment  sensitivity  with  respect  to 
A  and  I.  For  the  dM/dA  approximations  in  Figure  44,  the  second-order  method 
generally  performs  better  than  the  first-order  counterparts.  At  locations  where  the 
first-order  direct  approximation  is  grossly  inaccurate,  the  best  the  second-order 
approximation  can  do  is  match  the  more  accurate  of  the  two  first-order 
approximations,  in  this  case  the  adjoint  approximate  solution  using  modal 
approximation.  An  example  of  this  trend  can  be  seen  around  the  40Hz  and  80  Hz 
points  in  this  figure.  For  the  dM/dl  approximation  error  in  Figure  45,  neither  of 
the  first-order  approximations  is  grossly  inaccurate,  and  hence  the  second-order 
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approximation  performs  well.  The  sole  exception  of  this  is  around  the  40  Hz 
point,  where  the  Ritz  vector  approximation  resulted  in  errors  well  over  100%  and 
again  the  second-order  approximation  was  only  able  to  match  the  better  one  of  the 
two  first-order  approximations. 

For  pressure  sensitivity  approximations,  the  second-order  approximation  is 
undisputedly  the  most  accurate  and  most  consistent  performer.  As  shown  in 
Figure  46  for  the  errors  in  dP/dA  approximations  and  in  Figure  47  for  the  dP/dl 
approximations,  the  second-order  method  consistently  outperforms  both  first- 
order  methods  significantly  when  approximating  pressure  sensitivities  in  the 
given  range  of  frequency.  While  the  errors  for  one  or  both  the  first-order  methods 
exceed  100%  in  some  areas,  the  second-order  method  performs  far  better  with 
maximum  error  of  40%  and  20%  for  dP/dA  and  dP/dl  respectively,  and  with 
average  error  far  lower  than  that. 

4.3.1  Adjoint  based  on  Ritz  vectors,  direct  based  on  modes 

For  comparison,  the  large-scale  model  was  again  solved  with  almost  all 
the  same  parameters  as  previously,  with  one  change.  In  this  method  the  adjoint 
solution  was  approximated  using  a  modal  basis  of  50  modes,  while  the  direct 
solution  uses  1  Ritz  vector  for  approximate  basis. 

Unfortunately  for  this  method,  the  Ritz  vector  approximation  was  simply 
unable  to  approximate  the  solution,  resulting  in  very  large  errors  for  the  direct 
solution  approximation.  Since  the  second-order  approximation  relies  on  having 
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two  reasonably  accurate  first-order  approximations  to  further  increase 
accuracy,  naturally  it  would  at  best  be  only  as  accurate  as  the  better  one  of  the 
first-order  approximations  if  the  other  one  is  grossly  inaccurate.  In  extreme  cases, 
where  one  of  the  first-order  approximations  has  consistently  1000%  or  higher 
errors,  the  effect  of  this  becomes  significant  enough  that  it  causes  the  second- 
order  approximation  to  show  higher  inaccuracy  than  even  the  alternate  first-order 
approach.  An  example  of  this  case  is  illustrated  in  Figure  48  where  large  errors  in 
the  adjoint  approximations  result  in  also  large  errors  in  the  second-order  method. 
In  this  case,  the  first-order  direct  approximation  actually  performed  better  than  the 
second-order  approximation. 

4.4  Convergence  study  of  the  second-order  approximation  when  using  a 
modal  basis  first-order  approximation 

In  this  study,  one  of  the  methods  for  obtaining  one  of  the  first-order 
approximations  utilized  in  the  second-order  approximation  method  uses  natural 
mode  shapes  of  the  structural-acoustic  enclosure  model  to  serve  as  a  basis  matrix 
for  the  approximation.  Typically  in  first-order  approximations  using  modal  basis 
it  is  expected  that  as  more  modes  are  used  to  approximate  the  model,  the  more 
accurate  the  approximation  becomes.  The  cost  of  doing  this  is  the  increased 
computational  cost  involved  in  solving  a  larger  problem  as  the  reduced  basis 
matrix  becomes  larger  as  more  modes  are  used.  At  the  extreme  end  all  natural 
modes  will  be  used.  Then  the  approximation  is  then  no  longer  an  approximation, 


43 

but  becomes  the  exact  solution,  except  that  computationally  it  costs  much  more 
to  solve,  as  there  are  additional  computational  costs  involved  in  obtaining  the 
mode  shapes  for  the  basis  matrix  (as  well  as  multiplying  the  basis  matrices  against 
the  full-order  set  of  equations). 

A  question  then  arises.  Does  this  increase  in  accuracy  with  increased 
modes  used  in  the  first-order  method  translate  to  greater  accuracy  when  using  a 
second-order  approximation?  To  answer  this,  the  same  small  model  that  was  used 
in  section  4.2.1  was  used  again  with  the  same  experimental  parameters  except  the 
modes  used  to  approximate  the  adjoint  solution  were  incrementally  increased 
from  10  modes,  20  modes,  and  then  40  modes.  Figures  49  and  50  respectively 
show  errors  of  the  second-order  bending  moment  and  pressure  response 
approximations  as  the  modal  basis  is  increased.  For  both  illustrations  it  can  be 
seen  that  an  increase  in  modes  used  leads  to  a  reduction  in  error  and  hence  an 
increase  in  accuracy.  Furthermore,  an  increase  in  modes  used  also  gives  an 
increase  in  the  range  of  frequency  where  the  approximation  is  accurate.  For 
example  it  can  be  seen  that  in  order  to  get  reasonable  accuracy  past  180  Hz,  at 
least  40  modes  are  required.  This  is  expected  because  the  lower  number  of  modes 
prevents  mode  shapes  in  the  higher  frequency  ranges  to  be  included  in  the  basis 
matrix.  Care  must  then  be  taken  to  ensure  that  a  significant  number  of  mode 
shapes  are  included  for  the  frequency  range  of  interest,  but  not  so  many  that  it 
unnecessarily  affects  computational  time. 
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Figure  4:  Error  (%)  in  Bending  Moment  approximation  as  a  function  of  excitation 
frequency.  Direct  solution  approximated  by  1  Ritz  vector,  Adjoint  solution  approximated  by 
41  modes.  Small-scale  (275  DOF)  model,  g=0.15. 


Figure  5:  Error  (%)  in  Pressure  approximation  as  a  function  of  excitation  frequency.  Direct 
solution  approximated  by  1  Ritz  vector,  Adjoint  solution  approximated  by  41  modes.  Small- 
scale  (275  DOF)  model,  g=0.15. 
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Figure  6:  Error  (%)  in  dM/dl  approximation  as  a  function  of  excitation  frequency.  Direct 
solution  approximated  by  1  Ritz  vector,  Adjoint  solution  approximated  by  41  modes.  Small- 
scale  (275  DOF)  model,  g=0.15. 
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Figure  7:  Error  (%)  In  dM/dA  approximation  as  a  function  of  excitation  frequency.  Direct 
solution  approximated  by  1  Ritz  vector,  Adjoint  solution  approximated  by  41  modes.  Small- 
scale  (275  DOF)  model,  g=0.15. 
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Figure  8:  Error  (%)  in  dP/dl  approximation  as  a  function  of  excitation  frequency.  Direct 
solution  approximated  by  1  Ritz  vector,  Adjoint  solution  approximated  by  41  modes.  Small- 
scale  (275  DOF)  model,  g=0.15. 
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Figure  10:  Error  (%)  in  Bending  approximation  as  a  function  of  excitation  frequency.  Direct 
solution  approximated  by  51  modes,  Adjoint  solution  approximated  by  1  Ritz  Vector.  Small- 
scale  (275  DOF)  model,  g=0.15. 
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Figure  11:  Error  (%)  in  Pressure  approximation  as  a  function  of  excitation  frequency. 
Direct  solution  approximated  by  51  modes,  Adjoint  solution  approximated  by  1  Ritz  Vector. 
Small-scale  (275  DOF)  model,  g=0.15.  Adjoint  Error  is  off  the  chart. 
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Figure  12:  Pressure  response  approximation  as  a  function  of  excitation  frequency. 
Approximate  direct  and  adjoint  solutions  are  the  original  (reference)  adjoint  and  direct 
solution.  Small  model,  100%  change  in  design  variables  of  the  variable  element,  g=0.06 
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Figure  13:  Pressure  response  approximation  as  a  function  of  excitation  frequency,  zoomed  to 
area  of  interest.  Approximate  direct  and  adjoint  solutions  are  the  original  (reference)  adjoint 
and  direct  solution.  Small  model,  100%  change  in  design  variables  of  the  variable  element, 

g=0.06 


Figure  14:  Error  (%)  in  Pressure  approximation  as  a  function  of  excitation  frequency. 
Approximate  direct  and  adjoint  solutions  are  the  original  (reference)  adjoint  and  direct 
solution.  Small  model,  100%  change  in  design  variables  of  the  variable  element,  g=0.06 
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Figure  15:  Error  (%)  in  dM/dA  approximation  as  a  function  of  excitation  frequency. 
Approximate  direct  and  adjoint  solutions  are  the  original  (reference)  adjoint  and  direct 
solution.  Small  model,  100%  change  in  design  variables  of  the  variable  element,  g=0.06 


Figure  16:  Error  (%)  in  dP/dA  approximation  as  a  function  of  excitation  frequency. 
Approximate  direct  and  adjoint  solutions  are  the  original  (reference)  adjoint  and  direct 
solution.  Small  model,  100%  change  in  design  variables  of  the  variable  element,  g=0.06 
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Figure  17:  Pressure  response  as  a  function  of  excitation  frequency,  zoomed  to  area  of 
interest.  Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint 
solution  is  by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  100% 
change  in  design  variables  of  the  variable  element,  g=0.06 
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Figure  18:  Error  in  Pressure  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  100%  change  in 
design  variables  of  the  variable  element,  g=0.06 
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Figure  19:  Error  in  dM/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  100%  change  in 
the  design  variables  of  the  variable  element,  g=0.06 
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Figure  20:  Error  in  dM/dA  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  100%  change  in 
design  variables  of  the  variable  element,  g=0.06 


61 


Figure  21:  Error  in  dP/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  100%  change  in 
design  variables  of  the  variable  element,  g=0.06 
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Figure  22:  Error  in  dP/dA  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  100%  change  in 
design  variables  of  the  variable  element,  g=0.06 
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Figure  23:  Pressure  response  approximation  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  100%  change  in 
design  variables  of  the  variable  element,  g=0.06 
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Figure  24:  Pressure  response  approximation  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  50%  change  in 
design  variables  of  the  variable  element,  g=0.15 
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Figure  25:  Error  in  pressure  response  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  50%  change  in 
design  variables  of  the  variable  element,  g=0.15 
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Figure  26:  Error  in  bending  moment  response  approximations  as  a  function  of  excitation 
frequency.  Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint 
solution  is  by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  50% 
change  in  design  variables  of  the  variable  element,  g=0.15 
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Figure  27:  Pressure  response  approximation  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Zoomed  to  area  of  interest. 
Small  model,  50%  change  in  design  variables  of  the  variable  element,  g=0.15 
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Figure  28:  Bending  response  approximation  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Zoomed  to  area  of  interest. 
Small  model,  50%  change  in  design  variables  of  the  variable  element,  g=0.15 
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Figure  29:  Error  in  dM/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  50%  change  in 
design  variables  of  the  variable  element,  g=0.15 
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Figure  30:  Error  in  dM/dA  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  50%  change  in 
design  variables  of  the  variable  element,  g=0.15 
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Figure  31:  Error  in  dP/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  50%  change  in 
design  variables  of  the  variable  element,  g=0.15 
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Figure  32:  Error  in  dP/dA  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  is  the  original  (reference)  direct  solution;  the  adjoint  solution  is 
by  modal  approximation  using  33  left  and  right  eigenvectors.  Small  model,  50%  change  in 
design  variables  of  the  variable  element,  g=0.15 
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Figure  33:  Error  in  Pressure  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  10%  change  in  design  variable  I  of  the  variable  element,  g=0.06 
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Figure  34:  Error  in  Pressure  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  20%  change  in  design  variable  I  of  the  variable  element,  g=0.06 
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Figure  35:  Error  in  Pressure  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  30%  change  in  design  variable  I  of  the  variable  element,  g=0.06 
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Figure  36:  Error  in  Pressure  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  50%  change  in  design  variable  I  of  the  variable  element,  g=0.06 
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Figure  37:  Error  in  dP/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  10%  change  in  design  variable  I  of  the  variable  element,  g=0.06 
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Figure  38:  Error  in  dP/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  20%  change  in  design  variable  I  of  the  variable  element,  g=0.06 
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Figure  39:  Error  in  dP/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  30%  change  in  design  variable  I  of  the  variable  element,  g=0.06 
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Figure  40:  Error  in  dP/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  50%  change  In  design  variable  I  of  the  variable  element,  g=0.06 
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Figure  41:  Error  in  dP/dl  approximations  as  a  function  of  excitation  frequency. 
Approximate  direct  solution  and  adjoint  solution  is  by  first  order  Taylor  expansion.  Small 
model,  50%  change  in  design  variable  I  of  the  variable  element,  g=0.06 


Figure  42:  Pressure  response  approximation  as  a  function  of  excitation  frequency.  Adjoint 
solution  approximated  by  50  modes,  Direct  solution  approximated  by  1  Ritz  Vector.  Large- 
scale  (1543  DOF)  model.  Zoomed  to  area  of  interest. 
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Figure  44:  Error  (%)  in  dM/dA  approximations  as  a  function  of  excitation  frequency 
Adjoint  solution  approximated  by  50  modes,  Direct  solution  approximated  by  1  Ritz  Vector 
Large-scale  (1543  DOF)  model. 
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Figure  45:  Error  (%)  in  dM/dA  approximations  as  a  function  of  excitation  frequency 
Adjoint  solution  approximated  by  50  modes,  Direct  solution  approximated  by  1  Ritz  Vector 
Large-scale  (1543  DOF)  model. 
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Figure  46:  Error  (%)  in  dP/dA  approximations  as  a  function  of  excitation  frequency. 
Adjoint  solution  approximated  by  50  modes,  Direct  solution  approximated  by  1  Ritz  Vector. 
Large-scale  (1543  DOF)  model. 


Figure  47:  Error  (%)  in  dP/dA  approximations  as  a  function  of  excitation  frequency. 
Adjoint  solution  approximated  by  50  modes,  Direct  solution  approximated  by  1  Ritz  Vector, 
Large-scale  (1543  DOF)  model. 


Figure  48:  Error  (%)  in  dP/dl  approximations  as  a  function  of  excitation  frequency.  Adjoint 
solution  approximated  by  1  Ritz  vector,  Direct  solution  approximated  by  51  modes.  Large- 
scale  (1543  DOF)  model. 


%  Error 


89 


Figure  49:  Reduction  in  error  of  the  second  order  bending  moment  approximation  as  a 
function  of  excitation  frequency,  as  modal  basis  is  increased.  Direct  solution  approximated 
by  1  Ritz  vector,  Adjoint  solution  approximated  by  10,  20  and  40.  Small-scale  (275  DOF) 
model,  g=0.15 
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Figure  50:  Reduction  in  error  of  the  second  order  pressure  approximation  as  a  function  of 
excitation  frequency,  as  modal  basis  is  increased.  Direct  solution  approximated  by  1  Ritz 
vector,  Adjoint  solution  approximated  by  10,  20  and  40.  Small-scale  (275  DOF)  model, 
g=0.15 


Chapter  5:  Conclusions 
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Results  were  presented  for  various  approximation  studies  using  the  second 
order  approximation  method  for  a  coupled  structural-acoustic  finite  element 
enclosure  surrounded  by  flexible  walls.  Application  of  the  combined  direct- 
adjoint  solution  for  second  order  approximation  was  extended  from  structural 
dynamic  models  to  coupled  structural-acoustic  models,  and  a  variety  of  first  order 
approximation  methods  for  the  direct  and  adjoint  problems  were  explored. 
Bending  moment/stress  and  pressure  frequency  response  and  sensitivity  to  design 
variables,  critical  for  optimization  techniques,  were  studied,  and  improvements 
were  quantified  for  the  accuracy  of  the  combined  direct/adjoint  second  order 
approximation  when  compared  to  first  order  approximation  techniques.  A  system 
where  the  design  variables  evolve  relative  to  a  reference  system,  as  is  typical  in 
the  optimization  process,  was  examined.  Finally,  the  effects  of  increasing  the 
modal  basis  for  the  second  order  approximation  method  were  also  presented. 

It  was  shown  that  the  second  order  approximation  method  was  able  to 
significantly  improve  accuracy  over  first  order  methods  provided  that  two 
reasonably  accurate  first  order  adjoint  and  direct  approximations  are  available.  In 
cases  where  only  one  was  accurate  and  the  other  was  grossly  inaccurate,  the 
second  order  approximation  would  only  be  as  good  as  the  better  one  of  the  two 
first  order  approximations.  Some  first  order  approximations  can  be  as  simple  as 
using  the  actual  direct  and  adjoint  solution  for  the  reference  model,  which  costs 
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nothing  computationally  to  obtain  beyond  the  reference  analysis.  More 
sophisticated  approximations  such  as  using  Ritz  Vectors,  Taylor  expansion,  or 
modal  approximation  could  be  employed.  A  large  part  of  this  work  has  been  an 
effort  in  finding  a  reliable  combination  of  adjoint  and  direct  approximations  for 
use  in  the  second  order  method. 

The  combination  that  gave  the  most  consistently  accurate  results  was  the 
adjoint  solution  approximated  using  the  modal  approximation  with  the  direct 
solution  approximated  with  using  the  Ritz  vector  method.  Conversely,  the  worst 
performing  approximation  was  the  adjoint  solution  approximated  using  the  Ritz 
vector  method  with  the  direct  solution  using  modal  approximation.  This  was 
primarily  due  to  the  inability  of  one  single  Ritz  vector  (obtained  for  one 
configuration)  to  capture  the  behavior  of  approximate  pressure  responses  when 
the  natural  mode  frequencies  were  shifted  even  slightly  due  to  changes  in  design 
variables. 

Where  the  second  order  approximation  method  has  been  shown  to 
improve  accuracy  it  also  offers  significant  computational  time  savings  by 
reducing  large  sized  problems  in  the  order  of  hundreds  of  thousands  of  degrees  of 
freedom  to  a  much  smaller  problem  of  mere  tens  of  equations.  For  the  purposes  of 
optimization  techniques,  these  computational  time  savings  are  compounded  as  the 
optimizer  recalculates  the  problem  many  times  in  the  optimization  process. 


93 

Further  investigation  of  this  method  (particularly  the  ways  for  finding 
more  reliable  and  cost  efficient  first  order  approximations  to  be  used  in  the 
combined  direct/adjoint  second  order  approximation  method)  can  potentially 
further  increase  accuracy  while  maintaining  the  cost  and  time  saving  advantage. 

It  is  hoped  that  this  work  contributes  to  the  efficiency  of  solving  design 
oriented  structural-acoustic  problems  by  providing  an  additional  practical 
numerical  computing  tool,  as  well  as  point  to  additional  avenues  of  research  in  the 
field  of  order  reduction  of  large-scale  design  oriented  finite  element  modeling. 
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Appendix  A:  Derivation  of  Acoustic  Equations 
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Assuming  a  compressible  homogeneous  fluid,  the  homogeneous  initial  pressure 
field  (uniform,  constant), 


Po>Po>Uo  =  0 

For  small  disturbances: 

p  «  Pq ,  p  «  p0,u  «  1 


(A.l) 


(A.2) 


The  Thermodynamic  and  state  relations: 

For  an  isentropic  process  (adiabatic,  no  change  in  entropy), 


Po  +  P  _  Picai  =c 

(po  +  pY  C 0,o, J 


dp  = 


total 

dP iota!  o 


dp=>  p  = 


dp  total 

dPtotal  o 


p  =  cjp^- 


dpt 


total 


=  Ap,oJ'1  = 

P total  P total 


(A.5) 

where  c  is  the  speed  of  sound  in  the  undisturbed  fluid. 
Hence,  for  the  pressure  and  density  perturbations, 

P  =  cjp 


(A.3) 


(A.4) 


(A.  6) 


Equations  of  motion 

Continuity: 
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a 

(A.7) 

For  small  perturbations: 

P total  ~  A>  +  P 

(A.8) 

V  =V 

(A.9) 

Neglecting  second  order  terms, 

£  +  p.V-v-0 
a 

(A.  10) 

Momentum 

Euler  equations: 

DV  V pma] 

DT  ptota, 

(A.  11) 

For  small  perturbations: 

(a>  +  p)-j^  =  -v(a>  +p)= -vp 

(A.  12) 

{Po  +  p{^  +  (v  '  V)v  j  =  -Vp 

(A.  13) 

where 

/_  f  d  d  d  V 

(v  •  V)v  =  v,  —  +  v2  —  +  v3—  v 

^  ^f|  £2^2  ^2X3  ^ 

(A.  14) 

This  results  in  a  second-order  term  and  can  be  neglected.  The  perturbation 
momentum  (Euler)  equation  is: 
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eft  v 
po-^= -y? 

a 

Combining  the  continuity  and  momentum  equations, 

|+p,v-»  =  0 

a 

&  v 


(A.  15) 

(A.  16) 

(A- 17) 


Single-field  formulations 

Pressure  Formulation: 

Differentiating  momentum  equation  with  respect  to  time  (interchanging  the  order 
of  differentiations): 


d2p  f  dv 

a2  +  p°v\a 


=  o 


Substituting  from  the  continuity  equation, 

d2p 


a1 


•  +  /?0V' 


KPo 


=  0  =>  —t-  +  V2p  =  0 

a 2 


a  2 

and  express  °  E.  in  terms  of 

a 2 


dp  _  dp  dp  _  1  dp  d2  p  _  1  d2 p 

a  a  cK  a  a2  cj  a2 

Thus,  the  acoustic  pressure  equation  is: 


1  d2p 

2  -v2 


=  V2p 


(A- 18) 


(A.  19) 


(A.20) 


(A.21) 
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with  boundary  conditions  (from  the  Euler  equation): 

Vp  =  -p0  —  on  the  boundary,  when  velocity  perturbations  are  known  (moving 
0  dt 

boundary  walls  in  an  enclosure).  Alternatively,  the  pressure  disturbance 
/>(*, ,x2,x2,t)  can  be  defined  on  the  boundary. 

A  velocity  potential  formulation: 

Starting  with  the  perturbations  and  continuity  equations  again  we  apply  the  curl 
operator  to  the  Euler  equation: 


f  +  0 

dt 


Vx 


ft  yj 
Po^  =  -VP 
dt  ) 


Po 


d(V  xv) 


dt 


■  -V  x  Vp  =  0 


Thus, 

4^  =  0^(Vxv)=0 

dt 


(A.22) 

(A.23) 


(A.24) 


The  perturbation  velocity  field  is  irrotational. 

The  zero  constant  value  for  the  curl  of  the  velocity  perturbation  vector  is  due  to 
the  fact  that  at  time  zero,  before  any  excitation  is  applied,  there  are  no 
perturbations,  and  the  perturbation  velocities  are  zero. 


If  Vxv  =0 

Then  the  velocity  vector  can  be  expressed  as  a  gradient  of  some  potential 
function: 


v  = 
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The  scalar  function  <fi(Xl,x2,x3,t)  is  defined  as  the  velocity  potential. 
Substitution  into  the  continuity  equation, 

a 


+p0vv  =  o 


C/t 


f  d<f> 


0 


(A.25) 

(A.26) 


Recall  that  the  reference  density  is  constant  throughout  the  fluid  before 
perturbations.  Before  perturbations  begin,  the  perturbation  pressure  and  velocity 
are  zero.  Hence, 


d<j>  .  d<j) 

Po^  +  P  =  ^^>  P  =  -Po-r 

a  a 


(A.27) 


Returning  to  the  continuity  equation,  expressing  the  density  perturbation  in  terms 
of  pressure  perturbation,  and  then  substitute  the  previous  relationship  for  the 
pressure  in  terms  of  velocity  potential. 


dp  _  dp  dp  _  1 

~a~lfr~a~~p\J  a2 

substitution  into  the  continuity  equation  leads  to 


(A.28) 


c 


1  d2<f> 

2  a 2 


00 


vV 


(A.29) 


Boundary  conditions: 

When  velocity  on  the  boundary  is  known  (due  to  motion  of  the  boundary) 
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v  =  V^5  on  the  boundary,  or  alternatively,  if  pressure  disturbance  is  given  on 
the  boundary,  then 


=  (A.30) 

dt  P0 

Weak  forms  -  the  basis  for  Finite  Element  Analysis 
The  differential  equation, 


1  d2p 

C 2  a2 


=  V2p  ->  V2p- 


multiply  by  a  scalar  function  gp(x 
fluid, 


_L^  =  0  (A.31) 

a2 


l>x2,x3,t)  ar>d  integrate  over  the  volume  of  the 


jj^Pdn--Ljjj^dn  =  o 


“'cc  n 


(A.32) 


Review  of  some  vector  integral  relations: 

Gauss  Theorem  (for  a  vector  field  ^  in  a  volume  q  with  surface  (boundary)  x): 


JJJv  •  AdQ  =  jj A  ■  ndL  (A-33) 

fl  Z 

where  fj  is  a  unit  vector  normal  to  the  boundary  surface  and  pointing  out  of  the 
enclosed  volume  of  fluid. 


To  integrate  JJj SpV2p  ■  dCl  us*ng  Gauss  Theorem,  we  can  define  a  vector  field: 
n 


A  =  dp- Vp 


(A.  34) 


In  this  case 
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v  •  a  =  vsp  Vp  +  Sp-  y2p 

Thus,  from  Gauss’s  Theorem 

JJJ(v  •  l)iCl  =  JJJ(v^  Vp  +  Sp- V2p}n  =  Jj SpVp •  ndZ 

n  n  i 

Then 


JJJl (dp  ■  V2p}n  =  \\SpVp  ■  TidZ  -  JJjv^  •  VpdO. 
n  in 

Note  that  the  directional  derivative: 

dn 

substituting  into  the  integral  equation  above  we  get  the  final  weak  form: 

J/#§  ■ -  J1JV* ' — T  <«  =  0 

i  n  n 

Note  the  symmetry  in  this  form  if  we  use  the  same  spatial 
N  (xl}x2,x3) *°  approximate  the  two  fields  p  and  §p . 


If  we  use  any  approximation  of  the  form 


p(xnt)  =  {Nl(xi)  N2(x,)  ...  #„(*,.)« 


<h(0 

<h(0 

QniO 


(A.35) 


(A.36) 


(A.37) 


(A.  3  8) 

functions 


(A.39) 
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Vp(xnt)  = 


VSp(xi,t)  = 


dxx 

dx2 

dx3 

fM. 

dxx 

?ML 

dx2 

tW 

8x, 


few} 


tew} 


In  the  volume  integral,  then 


Sp{xi,t)  =  [Nx{xi)  N2(x,)  ...  iVn(x,.)U 


The  gradients: 

'w 


5qx{t) 

Sq2(t) 

Sqnif) 


Mis,} 


VSp  ■  Vp  =  {Sq}7 


a{N)  d{N]  d{N) 


dx. 


dx% 


£M. 

ckx 

m. 

dx2 

?¥L 

dx , 


fe(')} 


(A.40) 


(A.41) 


(A.42) 


(A.43) 


which  leads  to 


VSp  ■  Vp  = 

r,orpM<?Mr  |  d{N}d{N}T  1  dfr}d{N}T~ 
<3c,  8x,x  dn2  ck2  dc3 

=WM?(o} 

(A.44) 

where  [&]  is  a  symmetric  matrix. 

Similarly, 

M?«)} 

ot 

where  [m]  is  symmetric. 
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Appendix  B:  The  Triangular  Acoustic  Finite  Element 

For  a  linear  C°  acoustic  triangle,  with  nodes  1, 2,  and  3  in  counter-clockwise 
cyclic  order 

The  area  coordinates  are, 


*  =  7)  AT2(£,/7)  N3(c f,77)U 


*3j 


y  =  {N^,r])  N2(Z,tj) 
where 


M 


y2 


> 


N2(Z,r?)  =  Z 

N3(Z,tj)  =  tj 

The  Jacobian  of  the  transformation  is  then 


dx  dx 


J  = 


d£  dt) 
dy  dy 


Using  equations  B.l  through  B.3, 


(B.l) 


(B.2) 


(B.3) 


(B.4) 


+  x2 


(B.5) 


dy 

Tj=-y,+y2 
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jjx 2-X,  ^2-^1 

[*»-*.  y3-yi. 
and 

detJ  =  ((x2 -x,Xy3  ~yi)~{x 3  -x,Xy2  -y,))=  2^ 
where  ^  is  the  area  of  the  triangle. 

Two  area  integrals  have  to  be  evaluated  for  the  triangle: 

^8p(x,y)-p(x,y)dxdy 

A 

and 

jjvSp(x,y)  ■  Vp{x,y)dxdy 

A 

Let 


P  =  {NM>n)  N2{^tj)  N^,tj) 
and 


8p  =  {Nl{Z,r1)  N2{Z,n)  N3(4,tj) 


8px' 

Sp2  > 
Sp2 


(B.6) 

(B.7) 

(B.8) 

(B.9) 

(B.10) 

(B  .11) 

(B-12) 


Then 
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J fo{x,y)  - p{x,y)dxdy  = 

A 


'K 

Pi 

Jjfe  Sp2 

n2 

■{N,  N, 

N,} 

Pi 

A 

A. 

Pi. 

det  J  •  dtjdrj 


(B.13) 


(B.14) 


(B.15) 


Note  that  in  (B.13)  through  (B.15),  N  is  a  function  of  (£,  tj)  . 

The  following  integration  formula  for  integrals  over  triangles  using  area 
coordinates  can  be  used, 


fflfyd&ir 


M! 

(a  +  b  +  2)! 


2  A 


(B.  16) 


Integration  of  the  terms  of  the  [rri]  matrix  for  the  element  leads  to 


The  area  integral 
jjv  Sp(x,  y )  •  Vp(x,  y)dxdy 

A 

for  a  triangular  linear  pressure  element  is  obtained  as  follows, 
First  converting  from  x,y  coordinates  to  area  £,77  coordinates, 


Vp  = 


'dp 

dp  d%  dp  drj 

~K 

dJL 

dp 1 

lap] 

dx 

d 4  dx  drj  dx 

dx 

dx 

.  _  r-i 

dp 

>  —  « 

dp  d£  dp  drj 

K 

j  dp 

'  — «/  ' 

9p| 

[Qy. 

d£  dy  di 1  dy 

dy 

dy_ 

W 

Sp 

and  expressing  the  area  coordinates  as  functions  of  x  and  y, 

A Tl(x,y)  =  ^-{iy2  -y3)x+(*3  ~x2)y+(y3x2 -x3y2)} 

2  A 

N2(x,y)  =2-{{y,  ~yi)x  +  {x ,  -x3)y  +  {y,x3  -x,^)} 

2  A 

n3  (*>  y) = T7  {O', *1  )y + (y2*»  -  *2Ti  )} 

2v4 

where  from  B.8, 

2A  =  (x3  - x2\yx  - y2)~ (x,  - x2Xy3  " T2) 


To  express  Vp  = 
be  used 


dx 

dp 


dx  J 
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(B.18) 


(B.  19) 


(B.20) 


in  local  area  coordinates,  chain-rule  differentiation  has  to 


dp  _dp  d%  dp  drj_ 
dx  dl;  dx  drj  dx 

Since 

Z  =  N2(x,y) 


(B.21) 


(B.22) 


and 


Ill 


7]  =  N3(x,y) 


Vp  = 


[dp) 

[dp ] 

dx 

l 

’(vs-^i)  (yi-yi)\ 

54 

Up 

~  2  A 

{X2~Xl)_ 

dp\ 

[dtj  J 

Since 

/>  =  {l-<?-77  <f  T]}\ 

and 

Sp  =  {\-4-j]  4  ii 

then 

'dp' 


Pi 

Pi 

Pi 


SPx 

SPi\ 

Sp3 


dp 


[5rj 

ddp 


'-I  1  O' 

Pi' 

-1  0  1 

Pi > 

Pi. 

5f 

dSp 


[S77j 


'-1  1  0" 

-1  0  1 

►  = 

► 

Sp3 

(B.23) 


(B.24) 


(B.25) 


(B.26) 


(B.27) 


(B.28) 


Finally  substituting  into  equation  B.18,  it  is  shown  that  for  a  triangular  element 
JJV<^(x,y)  ■  VpO,  y)dxdy  = 

A 


T 

r-i  -n 

■■rJ) 

i 

0 

0 

1 r 

Sp, 

Sp2 

8Pi  J 

’(yj- 

(B.29) 

Thus,  for  a  triangular  element. 


lv3-^i)  fo“*J 


T_1  1  0 

d&i}~ 

Pi 

[-1  0  1 

Pi  ‘ 

Pi. 

112 


$P\ 

' Px 

\\y/Sp(x,  y )  •  Vp(x,  y)dxdy  =  < 

Sp2 

W- 

Pi ' 

A 

Sp, 

Pi. 

where 


(B.30) 


-1  -1 

1  0 

0  1 


rta-Ti) 

Ifa-yi) 


{xi~x3) 

{X2~Xl). 


~(y3-yi) 


Cvi  -  >^2  )T-  1  °* 

(x2 -1  0  1 


(B.31) 


It  is  also  noted  that 

(y3  -yj  +  ixi  -  )2  =  4 2  (B-32) 

Cvi  - y2 Y  +  (*2  ~xiY  =  Ln  (B-33) 

iy3  -TiX^i -jO+fa  -  xiXxi  -  xi)=^(.Li2  ~Li2  ~Ln)  (B-34) 


where 

Ltj  is  the  distance  between  vertex  i  and  vertex  j  on  the  triangle. 

For  a  triangular  acoustic  element  with  nodes  1,2,  and  3  counter-clockwise  at  its 
vertices,  the  vertex  coordinates  being  xl,ylx2,y2x3,y3  then 
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Ln  =iy3-yif +{xi~xJ 
Ln  =(y  i -yi)2 +{x2~xi)2 

l232  ={y3~y  if +(x2-h)2 

(B.35) 

and 

A =\((x2 -xifoi -yi)~ix3  ~xi\y2  -  j'i))  (B-37) 

Performing  the  matrix  multiplication  in  B.31  and  substituting  the  identities  in 
B.32-34,  the  stiffness  matrix  [k]  can  then  be  written  as 


M 


is2  2,2  +  iU2)  |(-i!.2+4S2-i,22) 

i,i2  ife2-!,,"-^2) 

T  ^ 

symm.  Ln 


(B.38) 


and  the  mass  matrix  as  given  in  B.17  is  then 


with  cl  being  the  speed  of  sound  of  the  fluid  concerned. 
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Appendix  C:  The  Structural-Acoustic  Coupling  Terms 


The  weak  form  of  the  acoustic  field  problem  is 

\\8p%-dZ-  jjjv dp ■  VpdQ  — \\\5p r/Q  =  0  (C.l) 

i  ™  n  coo  n  a 

On  the  boundary  between  the  acoustic  enclosure  and  structure  that  surrounds  it 
the  surface  integral  is 


f  \Sp%-dL  (C.2) 

i 

The  Euler  equation  for  small  perturbations  is 

Po~  =  ~^P  (C-3) 

dt 

Equation  C.3  makes  it  possible  to  relate  the  normal  derivative  of  pressure  on  the 
boundary  cpjdn  to  the  motion  of  that  boundary.  Multiplying  by  a  direction  vector 
normal  to  the  surface  (and  pointing  outwards), 


du  _  _  dp 

p0--n=-Vp-n=-— 
dt  dn 


thus 


(C.4) 


dn  0  dt 


(C.5) 


where  u  ■n  =  u„  is  the  velocity  component  normal  to  the  surface  in  the  outwards 
direction.  («  here  denotes  velocity,  not  displacement  as  in  the  structural  finite 
element  formulations). 
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For  a  triangular  acoustic  element  with  nodes  1-2-3  (counter  clockwise)  touching  a 
beam  element  A-B,  the  derivation  here  will  be  limited  to  the  case  where  the  nodes 
A,B  and  the  nodes  1,2  coincide,  respectively. 


The  surface  integral  in  this  2D  case  becomes  a  line  integral  from  A  to  B  (1  to  2) 

(c.6) 

For  a  2D  beam  element  in  local  coordinates  £,  7 , standard  FE  interpolation  for  the 
transverse  displacement  is  given  by 


»<f)  =  W.'(?)  K(4)  »l(S)  Kit)) 


w\ 

iex 

W2 

L6-, 


(C.7) 


where  for  j-  =  £  e  [0,l] : 

N^)=  1-3#2  +2#3 
N&)  =  {-  2{2+t3 
N&)  =  3<f-2£3 
n:(z)=-42+{3 


Along  the  same  line  (1,2)  the  acoustic  pressure  disturbance  is  linear 

N:{4ty-{*'  *6} (c9> 


Using 
dx  =  Ld£ 


(C.10) 
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we  have 

§Sp%dZ  =  ~P,‘\Sp~c =  SP'.'f  PJ-' 

j{^||}kfe)  n'As)  N’Ai)  iv.'fete- 

Note  that  there  is  a  change  of  sign  due  to  the  fact  that  the  transverse 
displacements  w  are  given  here  as  pointing  into  the  acoustic  enclosure,  while  the 
surface  integral  takes  as  positive  displacements  pointing  out. 

For  a  beam  element  oriented  with  respect  to  the  global  x-axis  by  an  angle  A ,  and 
including  both  axial  and  transverse  deformations,  transformation  from  local  to 
global  displacements  leads  to 


w, 


L9X 
L6 ' 


(C.ll) 


A  direction  vector  along  the  axis  of  the  beam 

fcosA 
sin  A 

where 

A  XB  ~XA 

cos  A  =  — - - 

L 

COsA  =  ^LlZl 
L 

and 


(C.  12) 


(C.13) 


L  =  tI(xb-xa)2  +(yB-yA)2 


(C.14) 
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In  global  coordinates  the  contributions  of  the  transverse  motion  w  to  motions 

Ux’Uy  is 


cA 

sA 


-sA 

cA 


cA 

-sA 


sA  ftO 

cAJrJ 


Then, 


w, 

-sA 

cA 

0 

0 

0 

O' 

L0, 

0 

0 

L 

0 

0 

0 

< 

w2 

0 

0 

0 

-sA 

cA 

0 

0 

0 

0 

0 

0 

Z 

Returning  to  equation  C.l  1,  we  now  have, 
jfo£-dZ  =  fo  5p2}-p,l)^^{n^)  N&)  N&)  N&)}d^ 


-  5  A 

cA 

0 

0 

0 

0" 

\i 

0 

0 

L 

0 

0 

0 

0 

0 

0 

-sA 

cA 

0 

«x  2 

0 

0 

0 

0 

0 

L 

»y2 

A. 

Expanding  the  term  from  C.l 7  above, 

KiS)  n!M)  N'At)}ds 


(C.l  5) 


(C.l  6) 


(C.l  7) 


we  arrive  at 


ifr(i-^i-3^2+2#3)  (i-^-2^2+#3)  (i-#X3£2-2<r3)  (1 -<?)(- £2 + <f 3 )' 

oJL  £(l-3<f2+2#3)  ^-2<f2+£3)  rf3£2-2#3)  <f(-£2+£3)  . 


Solving  the  integral  gives 
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And  the  matrix  multiplication 


-sA 

cA 

0 

0 

0 

O' 

[7/ 

720 

Yio 

3/ 

720 

l 

o 

< — 

1 

0 

0 

L 

0 

0 

0 

3/ 

L720 

Xo 

7/ 

720 

XoJ 

0 

0 

0 

~sA 

cA 

0 

0 

0 

0 

0 

0 

L 

gives 

f-  Is  A./  7  c  A/  L/  -3s  A/  3c  A/  -1/  ' 

/ 20  /20  /20  720  720  /30 

-3s  A/  3c  A/  1/  -Is A/  7  c A/  _L/ 

L  720  720  730  720  720  720 J 


Then,  the  contribution  to  the  surface  integral  of  a  beam  and  an  acoustic  element 
touching  it  is 


fy2}-PoL 


Is  A/ 

720 

2sA/ 

/20 


7  c A/  L/ 
720  /20 
3cA/  L/ 
720  /30 


-3s A/ 
720 

-Is  A/ 

720 


3cA/ 

720 

7  c  A/ 
720 


V 
$  , 

1 

N>< 

w\ 

o 

L/l  0. 

«*2 

Uy2 

/J 

1^2  J 

(C.18) 


Or  simply, 

\\8P%^{3pY.PMM-L  (C-19) 

I  ^ 

where  {qs  }is  the  vector  of  structural  degrees  of  freedom,  and  [Sp]  is  the  vector 


of  virtual  changes  in  the  pressure  degrees  of  freedom. 
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Combining  Structural-Acoustic  terms. 

Recalling  from  Appendix  B  that  for  the  area  integrals  for  a  single  acoustic 
element  are 


'8P\ 

T 

a' 

||V4?(x,y)  •  Vp(x,y)dxdy  =  • 

•W 

Pi 

A 

.Pi. 

Assembling  contributions  from  all  elements  we  get  the  integrals  for  the  complete 
acoustic  cavity 


\\V8p(x,  y)  •  Vp(x,  y)dxdy  =  {Sp}T  [Ka  ]{/>} 

A 

\  | ]8p{x,  y )  •  p(x,  y)dxdy  =  {< 5pJ  [Ma  ]{/>} 

CO0  A 

FE  discretization  of  the  acoustic  equations  in  weak  form  (from  C.l), 
\\Sp%dL-  ifrSp  ■  Vpdn  - -4  \\\Sp ^dQ  =  0 
now  leads  to 


{SpY  lMa  }{p}+[Kjp}- Po  [R„  U  }]  =  0  (C.20) 

Then 


lMa}{p}+[Ka]{p}-Po[RasU}]=  0 


(C.21) 
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For  the  structural  part  (the  beams,  assuming  viscous  damping  in  the  structure 
and  applied  external  forces)  we  have 

Ik  fe } + [c,  ]{?, } + [k,  \q, }  -  k  rwl=w  <c-22) 

The  pressure  contributions  to  the  structural  equations  are  obtained  by  integration 
of  pressures  over  the  length  of  each  beam,  pressures  providing  the  applied 
transverse  load  on  the  beam. 


Thus  the  final  coupled  structural-acoustics  equations  are 


Ms 

0  ' 

cf 

0  " 

H+ 

~K, 

_  PoRas 

Ma_ 

l^J 

0 

C..' 

[Pi 

0 

Ka  J| 

F 

0 


(C.23) 


Appendix  D:  The  Second  Order  Approximation 
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The  second-order  approximation  is  derived  from  a  process  comparable  to  central- 
differencing  second-order  derivative  approximations.  The  key  idea  is  to  eliminate 
first-order  error  terms  by  manipulating  alternative  expressions  for  the  desired 
response. 

Let  the  direct  and  adjoint  solutions  of  a  linear  matrix  equation  be 

MW  =  (M  (D1) 

M*W={c}  (D-2) 

where  response  y  can  be  calculated  from  either  the  direct  solution  {x}  or  the 
adjoint  solution  {r]}  by 

y  =  W  *  M  =  {c}  *  [A]'1  {&}  =  {?7}*  {*}  (D-3) 

An  error  {<&}  in  the  direct  solution  will  result  in  a  corresponding  error  in  the 
calculated  response 

Syx  =  {<?}*{&}  (D.4) 

and  likewise,  an  error  {Stj}  in  the  adjoint  solution  will  result  in  a  corresponding 
error 

dy2={Srj}*{b}  (D.5) 

If  there  are  no  errors  in  Eqs.  (D.l)  and  (D.2)  and  their  solution  {*}  and  { 77},  then 
the  {b}  and  {c}  vectors  in  Eqs.  (D.4)  and  (D.5)  can  be  replaced  by 
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Syi={tl}*[A{&c\  (D-6) 

(°-7) 

Adding  (D.6)  and  (D.7), 

+  Sy2  =  {77  +  St]}t[a]{x  +  &c}- {77}*  [^]{x}-  {Stj}*  [a]{&c]  (D.8) 

The  first  two  terms  are  the  error  in  the  expression  {77}*  when  the  exact  {*} 
and  { 77}  are  replaced  by  approximate  {*}  +  {<&}  and  {77}  +  {Sr/} .  Based  on 
equations  D.4-D.8, 

S{{c}  *  {*})  +  ^({77}  *  {b})  -  <?(M*  MM)  =  -{5rl} *  MM)  (D-9) 

Equation  D.9  is  the  error  in  the  expression 

{c}  *  {*}  +  {77}  *  {6}  -  {77}  *  [Ajx}  (D.  10) 

when  the  exact  direct  and  adjoint  solutions  are  replaced  by  approximate  ones.  The 
expression  in  D.10,  when  exact  direct  and  adjoint  solutions  are  used  together  is 
just  the  response  y  because  with  exact  solution  Eqs.  (D.l)  or  (D.2)  can  be 
substituted  into  Eq.  (D.10)  to  cancel  either  {77} *  {£>}-  {77}*  MM  or 

W*M-M*MM- 


